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Motivated by the recent "Cosmos Project" observation of dark-matter concentrations with no 
ordinary matter in the same place, we study the question of the existence of compact objects made 
of pure dark matter. We assume that the dark matter is neutralino, and compare its elastic and 
annihilation cross sections. We find that the two cross sections are of the same order of magni- 
tude. This result has a straightforward and important consequence that neutralinos comprising a 
compact object can not achieve thermal equilibrium. To substantiate our arguments, by solving 
Oppenheimer-Volkoff equation we constructed a model of the star made of pure neutralinos. We 
explicitly showed that the condition for the thermal equilibrium supported by the Fermi pressure 
is never fulfilled inside the star. This neutralino state can not be described by the Fermi-Dirac 
distribution. Thus, a stable neutralino star, which is supported by the Fermi pressure, can not 
exist. We also estimated that a stable star can not contain more than a few percents of neutralinos, 
most of the mass must be in the form of the standard model particles. 



I. INTRODUCTION 

Recent observational data [l], 0, Q indicate that the 
universe contains a significant fraction (22%) of dark 
matter whose origin is still unclear. A possible solu- 
tion to this problem comes from supersymmetric (SUSY) 
models in the form of neutralino [4J. Neutralino is the 
lightest supersymmetric partner in SUSY, with the mass 
of about lOOGeV, and is stable. It interacts with the 
gravitational and weak interactions only which indicates 
that it is " dark" . Weak interactions and neutralino mass 
are sufficient to satisfy the relic density needed to explain 
the observed portion of the dark matter in the universe. 

In a remarkable recent study a detailed distribu- 
tion of dark matter function of the redshift in a 
part of our universe was obtained. These observations 
indicate that dark matter plays a role of a scaffolding 
upon which ordinary matter builds structures. However, 
the observations show that large pockets with only dark 
matter (and no ordinary matter) also exist. Then a natu- 
ral question arises whether compact objects like planets, 
stars or maybe even larger structures made exclusively 
of dark matter can exist. 

Recent study discussed the existence of stars powered 
by annihilation of neutralinos However, these ob- 
jects contain mostly baryons and only several percents of 
dark matter. Compact objects made of pure dark mat- 
ter were discussed in 0, H,M G3 ■ There it was assumed 
that the equilibrium of these stars was supported by the 
Fermi pressure. It was argued that the size of such a 
dark matter star is less than Im with the mass less than 
10~ 4 M Q , if neutralino mass is around lOOGeV Q. Since 
the density inside a neutralino star is very high, neutrali- 
nos collide frequently. These earlier studies used the esti- 
mate that identical non-interacting fcrmions have elastic 
cross section of about ir\ 2 , where A is de Broglie thermal 
wavelength With this assumption, non-relativistic 

neutralinos have the clastic cross section which is much 



larger than the weak interaction cross section. If this 
claim is indeed correct, than a thermal neutralino star 
could exist, since neutralino annihilation can be neglected 
in this case. However, it is unlikely that the neutralino 
clastic cross section is as large as nX 2 . For example, the 
large elastic cross section would directly contradict the 
recent "bullet cluster" observations It is therefore 

very important to explicitly check this fact. 

In order to achieve a thermal distribution, the lifetime 
of a particle has to be much longer than the characteris- 
tic interaction time. Stable neutralinos could in principle 
achieve thermal state in cosmological context. However, 
at high densities annihilation can not be neglected. A 
necessary condition for the thermal equilibrium among 
neutralinos (at densities smaller than (lOOGeV) 4 ) is that 
the ncutralino-ncutralino clastic cross section is much 
larger than the neutralino-ncutralino annihilation cross 
section. 



II. NEUTRALINO ELASTIC AND 
ANNIHILATION CROSS SECTIONS 

In this section, we will explicitly calculate and compare 
the neutralino-neutralino elastic and annihilation cross 
sections. We will show that they are comparable, i.e. of 
the same order of magnitude. This implies that neutrali- 
nos can not be in a thermal equilibrium in a compact 
object. After several encounters they will annihilate and 
be converted into the standard model particles. The ther- 
mal state of neutralions can not be reached. Thus, the 
neutralinos can not be described by the the Fermi-Dirac 
distribution. This then indicates that compact objects 
made only of neutralinos and supported by the Fermi 
pressure can not exist. 

In SUSY, the neutralinos Xi are the mixed states of the 
neutral higgsinos, '5^ and ^ 2 H , and the neutral gaug- 
inos, binos (Xb) and winos (X\) [l3[. The mixing and 
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mass matrices are given by 
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FIG. 1: The possible vertices of neutralino-neutralino inter- 
action. 
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where Zjq is the matrix that diagonalizes the mass ma- 
trix. Parameters Mi, M 2 and \i are mass parameters 
in the potential, v\ and v 2 are the vacuum expectation 
values of the two higgsions, Sw and Cw are sinus and cos- 
inus of the weak angle, while e is the electron charge. The 
matrix Zn also defines the mixing between the states 
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Here, M v o < M v o < M v o < M v o. 
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Neutralino x? is presumably the lightest supersymmct- 
ric particle, and can not decay into other supersymmct- 
ric particles. Because of the R-parity conservation, it 
can not decay into the standard model particles either. 
Therefore, neutralino is stable. The only way to destroy 
ncutralinos is annihilation. 

The cross section for the two-particle final state is 
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In our case, both particle A and particle B are neutrali- 
nos. In an elastic collision, particles 1 and 2 are neu- 
tralinos. In an annihilation process, particles 1 and 2 are 
the standard model particles. Pi is a momentum of the 
particle 1 . va and vb are velocities of particles A and B 
respectively. Ea and Eb are energies of particles A and 
B respectively. M is the scattering amplitude. If the two 
incoming ncutralinos have non-relativistic velocities, say 
v ~ lOOkm/s, then the momentum Pi of the light stan- 
dard model particle (except maybe for the top quark) in 
an annihilation process is about v/c ~ 1000 times larger 
than the momentum Pi of the neutralino in an elastic 
collision. Thus, only for highly relativistic incoming neu- 
tralinos clastic collisions will have Pi of the same order of 
magnitude as an annihilation process. Since for various 
different channels all of the parameters are fixed, except 



for the amplitude, in what follows we compare the am- 
plitudes only. 

We now look for the possible neutralino-neutralino 
elastic interactions. The possible neutralino-neutralino- 
something vertices can be found in [l3| and are shown in 
Fig. [T] The values of the vertices are 
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Here, mz 
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(vf + v 2 ) and m w 



v 2 ) ' 5 . The masses of the Z and W bosons are mz and 
mw respectively. From these vertices, the possible Feyn- 
man diagrams for clastic neutralino-neutralino interac- 
tions are shown in Fig. [21 Consider now non-relativistic 
neutralino-neutralino elastic collision. The amplitude in 
the s-channel is approximately 
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To simplify these equations, omit writing the projectors 
Pl and Pr, since they will not affect the order or mag- 
nitude estimate. Here, s ~ 2M 2 , while m is the mass of 

Hi and A° k . M^f H0 ^ is the amplitude if H% or Al are 

the mediators, while M^' z0 is the amplitude if Z® is the 
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FIG. 2: The neutralino-neutralino elastic collisions: The 
dashed line can be //J?, and From the top down 

are s-, t- and u-channel respectively. 



mediator. The t- and u-channcl clastic amplitudes are 
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In this case 



0. If masses of 



H®, and are comparable to M x o (and the energies 
are not tuned to be in the resonant channel) , then all the 
channels are of the same order of magnitude. 

We now compare the s-channel elastic scattering am- 
plitude with the s-channcl annihilation amplitude (shown 
in Fig. [3]). For the simple order of magnitude estimate we 
ignored the production of Z and W bosons. We note that 
this a conservative approximation since these channels 
will only enhance conversion of neutralinos into standard 
model particles. In order to find the annihilation am- 
plitude we need to calculate the three vertices in Fig. |U 
The values are 
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Parameters D, E, Yf and j have different values for dif- 
ferent fermions. For it-quark we have D = — 1, E = g, 
Yf = \/2mf/v2 and j — 2. For d-quark we have D = 1, 
E = §, Yf = \f2mfjv\ and j = 1. For e we have D = 1, 
E = 2, Yf = \/2m f/v\ and j = 1. For v we have D = — 1, 



0, F/ = and j = 2. The amplitudes are 
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where f c is antiparticle of /. Again we omit writing Pj, 
and Pr. Eq. p5)l gives the amplitude if the mediator is 
while Eq. (fl~6|) gives the amplitude if the mediators 
The amplitude in Eq. (fT5]l is of the same 
order of magnitude as the amplitude for the elastic scat- 
tering in Eq. ||5J). The amplitude in Eq. (|16p is roughly 
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FIG. 3: The neutralino-neutralino s-channel annihilation: 
The dashed line represents Zz, H®, and A°. 
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FIG. 4: The s-channel annihilation vertices. 



(rrif/Mz) times the amplitude for the elastic scattering 
in Eq. @ . If m/ is of the same order of magnitude as 
Mzi we can conclude that the clastic and annihilation 
scattering amplitude are also of the same order of mag- 
nitude. The only case when this may not happen is when 
the neutralino is lighter than the top quark and its ini- 
tial energy is very low. In that case top quark can not 
be in the final state for kinematic reasons. However, as 
mentioned earlier, the cross section also includes the mo- 
mentum of the emitted particle. Since the next particle 
in the mass hierarchy is b-quark, which is much lighter 
than neutralino, the final state particle will have a large 
momentum. Then the annihilation cross section will be 
enhanced by a large factor as explained after Eq. ([?]). 
Also, there arc more branches in the annihilation process 
than in an elastic collision since the fermion / in Fig. [3] 
stands for several species of particles. This is more than 
enough to overcome the mass suppression factor if the 
neutralino is lighter than the top quark. Finally, for our 
argument to work, it is enough that there exists only one 
annihilation channel which is comparable to the elastic 
cross section. As we will show, t- and u- annihilation 
channels have no mass suppression factors regardless of 
the sign of the neutralino-top quark mass difference. 

We now compare the t- and u-channcls elastic scatter- 
ing amplitude with the corresponding annihilation am- 
plitude (shown in Fig. [5]). In this case, the mediators are 
sfermions. The vertices needed for these diagrams are 
shown in Fig. [5] The values are 
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where Af and Bf are constants, which are different for 
different species [13]. The amplitude is 
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FIG. 5: The neutralino-neutralino t- and u-channel annihila- 
tion: The dashed line represents a sfermion. 
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FIG. 6: The vertices of fermion-neutralino interaction. 



Again we omitted writing Pj, and Pr. Here, v 
square of the momentum transfer, while is the mass 

of sfermion /. If rrij is not much different from the mass 
of Z°, H®, and A^, then the t- and u-channel annihilation 
amplitudes are of the same order as the t- and u-channel 
elastic amplitude in (JTUJ) . 

This analysis indicates that the neutralino annihila- 
tion amplitude is at least of the same order as the neu- 
tralino elastic amplitude. The annihilation cross section 
can be much larger than elastic cross section, when one 
includes the momentum difference between the light stan- 
dard model quarks and heavier neutralinos in the final 
state, as explained after Eq. ((4|. 



III. MODEL OF THE STAR MADE OF PURE 
NEUTRALINOS 
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FIG. 7: The black dots are particles moving with velocity 
v(E). A small wall of area A is embedded inside the particle 
distribution. 



transfer of particles that hit the area A. In a three- 
dimensional space, about 1/6 of the particles move to- 
ward one side of the wall and bounce back. The momen- 
tum transfer, V, in the time interval At is 



AT(E) = 



n(E)AEv 2 At 



The pressure is 



P = 



AT(E) 
AAt 
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nEv 2 
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(22) 



This is a common result in statistical mechanics. Note 
that the crucial assumption in Eq. (|2ip is that every par- 
ticle that arrives at the wall A is bounced back. In other 
words, the efficiency of the interaction between the wall 
and particles is 100%. If the wall bounces back only j(E) 
of the incoming particles (i.e. 1 — ~f(E) of the particles 
pass through the wall) Eq. (f2"Tj) must be rewritten as 



AT(E) = 1 {E) 



n{E)AEv 2 At 



(23) 



In this case, the pressure must be rewritten as 



In this section we quantify our arguments about the 
non-existence of the star made of pure neutralinos by 
treating the Pauli exclusion force as a real interaction 
between the fermions (different than the weak interac- 
tion) . 

First we will show that the mean free path for the 
interaction caused by the Pauli exclusion force must be 
shorter than the distance between fermions if the Fermi 
pressure is to provide stability. To show this, we consider 
a general setup like in Fig. [7] A wall with an area A 
is embedded in a group of particles (either fermions or 
bosons). The number density of particles with energy 
E is n(E). The particles move with velocity v(E). The 
total particle number density, N, and energy density, p, 
are 



N = J n(E)dE 



(19) 



p = J n(E)EdE (20) 
The pressure can be calculated from the momentum 



AV{E) 
AAt 



dE = 



j(E)nEv 2 



(24) 



The factor ^(E) then measures the departure from the 
thermal pressure {"l(E) = 1 means the pressure is equal 
to the thermal pressure). 

Now consider the wall made of particles themselves. 
From the number density, the size of a particle is I = 
n -i/3 rp ne cruc j a j factor 7 (£7) can be estimated from 
the ratio between the interaction cross section between 
the particles, <j(E), and the area I 2 : 



I 2 
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Here, A = \/{na) = l 3 /<r is the mean free path of the 
particles. If A < I, i.e. the mean free path is shorter than 
the distance between the particles, then we need to set 
j(E) — 1. However, if the mean free path is larger than 
the distance between the particles, j(E) < 1. In that 
case, the pressure is smaller than the thermal pressure. 
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In the extreme case when A I the particles practically 
do not interact with each other, i.e. j(E) w 0, and the 
pressure can not be built up at all. These arguments 
are true both for bosons and fcrmions. For fermions, 
the Pauli exclusion force builds up the Fermi pressure 
between the fermions. However, for this to happen, the 
mean free path for fcrmions must be smaller than (or at 
most equal to) the distance between the fermions. We 
will use this result in what follows. 

We now model the structure of the star made of pure 
dark matter. The relativistic mean field approximation 
[lH implies that the density, pressure and number density 
of neutralinos can be written as (for details see fioj ) 
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FIG. 8: The dependence of the ratio between the mean free 
paths for the Pauli exclusion force and annihilations \ p /\ a on 
the Fermi momentum kF- The condition (|28fl breaks down for 
high kF (high densities). 
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Here, kp is the Fermi momentum. Also 
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m* = m x + gnH + g^h 

The numerical values of the relevant parameters are 
m x = 1.01679763 x l0 e MeV, m H = 7.5336438 x 
10 5 MeV, m h = 1.18760345 x 10 5 MeV, m z = 9.1187 x 
10 4 MeP, g H = 3.131 x 10" 1 , g h = -2.705 x 10~ 2 , and 
g z = 1.43 x 10~ 3 . 

The ncutralino annihilation cross section is a a ~ 1 x 
10~ 10 MeV~ 2 . The mean free path for annihilations is 
A a = l/(Ncr a ). As we showed above, the mean free path 
of the Pauli exclusion force is about the distance between 
fermions, i.e. \ p = N^ 1 / 3 . To keep neutralinos in ther- 
mal equilibrium, we need 



^ = N 2 ' 3 a a « 1 



(28) 



can not 



From Fig. [51 we see that the condition 
hold at high kp (i.e. high densities). 

We now write the Oppenheimer-Volkoff equation [l 



dP 

~dV 

dM 
dr 



(p + P){M + 47rPr 3 ) 



r 2 (l - 2M/r) 



(29) 
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Combining Eqs. (|29|) with the equation of state given by 
Eqs. (ptJ)) . we can find the structure of the neutralino 
star. We need to setup the boundary condition for /cp, 
so we choose kp = 0.5m x at the center of the compact 
object. We then numerically integrate Eqs. (|29[) . The 
radius of the star is the distance from the center at which 
the pressure drops to zero (in this case the radius of the 
star is 1.2cm). This gives us the structure of the star, i.e. 
energy density, pressure and neutralino number density 
as functions of the distance from the center of the star. 

Fig. [S] shows the ratio of X p /X a vs. radius. Clearly, 
the ratio is about 1 throughout the star. The ratio drops 
sharply only at the surface of the star. This indicates 
that the Pauli exclusion force is not sufficient to overcome 
annihilation. The neutralinos decay into the standard 
model particles before the star becomes thermal. 

One may roughly estimate the ratio between the stan- 
dard model particles and neutralinos which is needed to 
stabilize the star. This ratio will depend on the mean 
free path for neutralino annihilation and that of the stan- 
dard model particles. For a condition (|28|) to hold, the 
mean free path of the standard model particles can not 
be greater than several percents of the mean free path 
for neutralino annihilation. This implies that the neu- 
tralino content in the star can not be greater than a few 
percents. 

Finally, we comment on the interplay between gravity 
and annihilations. Gravity brings neutralinos together, 
but annihilations put an upper limit to the density of 
(pure) dark matter clumps. Even without detailed nu- 
merical simulations we can perform a simple order of 
magnitude estimate. Consider for example a neutralino 
clump that was presumably created when the galaxy was 
formed, and survived till now. The density of the clump 
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FIG. 9: The dependence of the ratio between the mean free 
paths for the Pauli exclusion force and annihilations \ p /\ a on 
the distance from the center of the compact object made of 
pure dark matter. Throughout the star the condition (|28|) is 
violated. The ratio \ v /\ a drops sharply only near the surface 
of the star. 



can be calculated from the condition 
AT<^ = 



<r„Nv^ 



(30) 



where AT is the lifetime of a galaxy — about 10 Gyrs, 
while the velocity of dark matter is assumed to be 100 
km/s. N is the total number density of neutralinos 
and the ncutralino annihilation cross section is a a ~ 
1 x lO~ 10 AIeV~ 2 . From this condition we find that 
the number density of the clump must be lower than 
N ~ 6.2 x 10" 35 GeV 3 « 8.2 x 10 12 uT 3 , which for a 
lOOGeV neutralino gives the mass density of 10 _13 g/cm 3 . 
For more general constraints, one would need to perform 
detailed numerical simulations of the gravitational col- 
lapse with annihilations included. 



IV. CONCLUSIONS 

In conclusions, we studied here whether pure dark 
mater can make compact objects like stars. The stud 
was motivated by the remarkable recent observations 
that there are regions in our universe with pure dark 
matter distribution without the visible presence of the 
ordinary matter. We explicitly calculated and compared 
the neutralino elastic cross section with the annihilation 
cross section and found them to be comparable. This 
means that, in each encounter, neutralinos have about 
the same chance to annihilate into the standard model 
particle or to be scattered. If they are scattered, they ex- 
change energy and momentum and after many repeated 
encounters will reach the thermal equilibrium. However, 
annihilation will prevent this scenario from happening. 
Once ncutralions get converted into the standard model 



particles they can not come back. This situation can not 
be described by the Fermi-Dirac equation of state. This 
implies that a stable neutralino star supported by the 
Fermi pressure can not exists. 

To substantiate our arguments, by solving 
Oppcnhcimer-Volkoff equation we constructed a model 
of the star made of pure neutralinos. We explicitly 
showed that the condition for the thermal equilibrium 
supported by the Fermi pressure is never fulfilled inside 
the star. We also estimated that a stable star can not 
contain more than a few percents of neutralinos, most 
of the mass must be in the form of the standard model 
particles. This may have implications for the question 
why we are made of the standard model particles if 
most of the matter in the universe comes in the form of 
dark matter. For intelligent observers like us to evolve, 
structures like stars and planets are (perhaps) necessary. 

One of the possibilities to have neutralinos in thermal 
equilibrium is an environment with enormous densities, 
of the order of (lOOGeV) 4 . If such an environment is 
thcrmalized, then its temperature is of the order of the 
rest mass of the neutralino. In this case neutralinos can 
be produced thermally. This can only happen in early 
universe or in the central core of a very dense star. In the 
case of the central core of a star, this re gion must be very 
tiny, with the radius smaller than 1cm [14| . Again, such 
an object can not be supported by the Fermi pressure. 
The annihilation process would be so quick that the star 
could not exist long enough to be observed. 

The other possibility to have neutralinos in thermal 
equilibrium is an environment with very low densities. 
In an environment where the neutralino density is so low 
that neutralinos barely collide with each other, their ef- 
fective lifetime is much longer than the weak-scale char- 
acteristic collision time. In such conditions, gravitational 
cooling may bring neutralinos to equilibrium. This is the 
mechanism behind dark matter halo formation. How- 
ever, these are rather low density dark matter distribu- 
tions and can not be called compact objects. At higher 
densities gravitational scattering can not compete with 
weak interactions. Therefore, gravitational cooling may 
bring neutralinos together but can not play an important 
role in the subsequent evolution of the compact object. 

We assumed here that the dark matter is comprised of 
neutralinos. However, any particle with the mass around 
the weak energy scale and weak scale interaction cross 
section (such is neutralino) automatically has a relic den- 
sity which satisfies observational constraints. Thus, we 
expect that the same conclusions will hold for more gen- 
eral models. Therefore, despite the fact that regions in 
our universe with the pure dark matter distribution were 
observed, compact objects like planets, stars and larger 
systems like galaxies made of dark matter are unlike to 
exist. 

We also note that we used the dark matter models in 
their most standard form. In some exotic models, dark 
matter particles with an extra self-interaction [rH ufL [l9l. 
I20I I21I [22l . I23I , [23 | were considered to solve the problems 
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with subhalocs [25|,|26||, cuspy cores [27 J2& [29|, |30| and 
recent apparent PAMELA [3l|, ATIC[Uand FERMI 
[33l ] observations. In the context of these models, this 
extra self-interaction has to be properly included. 
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